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EFFECTS OF SUBCOOLING AND GRAVITY LEVEL ON BOILING 
IN THE DISCRETE BUBBLE REGION 
by Thomas H. Cochran and John C. Aydelott 
Lewis Research Center 

SUMMARY 

The effects of changes in subcooling and gravity level are investigated for water boil- 
ing at low heat flux from a flat horizontal surface. The zero gravity data were obtained 

by allowing the experimental package to free fall in a 100- foot drop tower, which per- 
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mitted the attainment of less than 10 times Earth gravity. An analysis of the dynamics 
of a generated bubble resulted in a new approach to the formulation of the buoyancy force 
and the definition of a force associated with the unbalanced pressure across the top sur- 
face of the bubble. The force associated with liquid inertia was obtained by the use of 
Newton's second law of motion. Application of the analysis to experimental data indicated 
that this newly defined pressure force was of major importance in bringing about bubble 
separation. Additional investigation showed that nucleate boiling was independent of grav- 
ity at high subcooling. 


INTRODUCTION 

The coming of space exploration has prompted investigations of the effect of gravity 
level on the transfer of heat to fluids. A particular area of interest has been the storage 
of cryogenic fluids in closed containers over long periods of time, as would occur on a 
coasting space flight. Under such conditions, the fluid may be subject to zero gravity 
subcooled nucleate boiling because of a relatively low heat flux from solar radiation. 

In recent years, a great deal of experimental work with nucleate boiling of saturated 
fluids in various gravity fields has been conducted (refs. 1 to 17). These investigations 
have shown that, in the low heat flux region, the boiling mechanism is dependent on the 
dynamics of vapor bubbles formed on the heated surface. Experimental studies of the ef- 
fect of subcooling on nucleate boiling in normal gravity by the University of Denver 
(refs. 18 to 20) also indicate this dependence. Analytical works , such as references 21 



to 24, propose the importance of bubble dynamics. 

The present work is concerned with the effect of subcooling and gravity level on the 
dynamics of steam bubbles formed on a flat horizontal surface. High-speed motion pic- 
tures were taken at subcoolings of approximately 5°, 10°, 15°, 25°, 35°, and 40° F in 
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normal gravity and in an environment of less than 10 times normal gravity. The heat 
flux selected was kept constant for all tests and set at a level in the range defined in nor- 
mal gravity as the discrete bubble region (ref. 13). In such a range in normal gravity, 
the bubbles generally grow undisturbed on the heated surface, and the heat-transfer mech- 
anism is highly dependent on the disturbance in the liquid caused by the growth and depar- 
ture of the bubbles. 

Analysis of the photographic data included measurements of growth characteristics of 
the bubbles and calculation of the forces acting on the bubbles during their growth. The 
formulation for the forces included a new treatment of the buoyancy force and the defini- 
tion of a force concerned with the unbalanced pressure across the top surface of the bub- 
ble. An original approach to the force associated with liquid inertia resulted from the 
use of Newton’s second law of motion. From force histories, interpretations regarding 
the dynamic effects are made and discussed. 

Motion picture supplement C-246, which shows the experimental procedure and sub- 
cooled boiling of water at the various test conditions, has been prepared and is available 
on loan. A request card and a description of the film are included at the back of this re- 
port. 

SYMBOLS 

A area, ft^ 

C Dr drag coefficient 

d diameter, ft 

Fg buoyancy force , lb force 

F Dr drag force, lb force 

Fgy dynamic force, lb force 

Fp pressure force, lb force 

F^ total removal force, lb force 

XX 

F c total surface tension force, lb force 
b 

F vertical component of surface tension force, lb force 

s y / 2 

g acceleration due to gravity, ft/sec 
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g c gravitational constant, (lb mass/lb force)(ft/sec ) 

2 

gQ standard acceleration of gravity on Earth, ft/sec 

M mass, lb mass 

P pressure, lbforce/sqft 

R radius, ft 

r general radius, ft 

T radius of curvature of top surface, ft 
t time, sec 

V bubble total volume, ft^ 

3 

bubble volume directly over base, ft 

Y bubble center of mass, ft 

y distance above heater surface, ft 

v velocity, ft/sec 

q 

p density, lb mass/ft 

2 

p dynamic viscosity, (lb force)(sec)/ft 

a surface tension, lb force/ft 

<p contact angle 

Subscripts: 

b base 

i inside bubble 

SL liquid 

max maximum 

0 outside bubble 

S heater surface 

sat saturated vapor 

v vapor 

1 location one 

2 location two 



ANALYSIS 


Bubble Model 

In order to calculate the volume and center of mass of the generated vapor masses, 
a bubble model is assumed. The only restrictions placed on the model are that it is sym- 
metric with respect to the y-axis, as shown in figure 1, and that the surface of the bubble 
directly over the base is spherical. The general nature of this model is prompted by 
Warner (ref. 18), who shows that, for subcooled conditions, specification of a particular 
bubble geometry, such as a sphere, is in error. 

Simple geometric formulas cannot be used to calculate bubble total volume because 
of the general nature of the model. Consequently, volume is determined by an integrative 
method in which a bubble is divided into segments, the volumes of which can be approxi- 
mated if they are assumed to be circular disks. The sums of the volumes of the disks 
produce the total volume. The volume directly over the bubble base is obtained by as- 
suming that the volume is a right circular cylinder with a segment of a sphere as a cap. 

The center of mass of a bubble is selected as the characteristic dimension to de- 
scribe its motion and is determined from the location of the plane parallel to the heater 
surface, which divides the bubble in half with respect to its total volume. 



Bubble Forces 

An analytical investigation of the dy- 
namics of a steam bubble on a heated sur- 
face results in the identification of buoy- 
ancy, surface tension, pressure, drag, 
and dynamic forces. The force associ- 
ated with the vapor weight is neglected 
because the experimental conditions are 
such that it is very small. However, for 
all fluids near the critical thermodynamic 
state and for some fluids that do not have 
a large difference between liquid and va- 
por density, such as hydrogen, this force 
is large enough to be included. 
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Buoyancy Force 


Buoyancy on an object submerged in a liquid is caused by the difference between the 
external hydrostatic pressure force on its top surface and the external hydrostatic pres- 
sure force on its bottom surface. Therefore, for a bubble attached to a heated surface, 
the volume of the bubble directly over the base does not generate a buoyant force because 
of the lack of liquid beneath the base. The buoyancy force, as derived in appendix A, is 
defined as 

f b - < v - V bW f- < l > 

g C 

and it acts in a direction opposite to the acceleration due to gravity. 


Surface Tension Force 


A surface tension force is generated at the boundary of a liquid and some other sub- 
stance, such as a vapor or a solid. Therefore, such a force exists at the base of a bub- 
ble attached to the heater surface at the boundary of the liquid, vapor, and solid surface. 
The direction of this force is perpendicular to the boundary and in the plane of the liquid- 
vapor interface. Surface tension is defined as the ratio of the surface force to the length 
along which the force acts, which, for an attached bubble, is 


a S “ 



( 2 ) 


The horizontal component of the force is cancelled out around the circumference of the 
bubble base and only the vertical component remains, such that 


F s = CT q 7rd b sin <p (3) 

y 

where the contact angle cp is the angle at the base between the liquid- vapor interface and 
the heater surface, as shown in figure 1. This force retards the movement of a bubble 
from the surface. 


Pressure Force 

The net force due to the uniform internal pressure on a bubble surface is zero for 
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a bubble surrounded by liquid. For a bubble attached to a surface, however, the net in- 
ternal pressure force on the spherical surface area directly over the base is unbalanced. 
Since the pressure is greater inside the bubble than outside, the force acts to remove the 
bubble from the surface. The detailed derivation of this force in appendix B yields 



Drag Force 

The movement of a mass through a fluid creates a retardation force, which is termed 
drag. Only an approximation can be made for this force because of a lack of knowledge 
regarding the motion of the fluid about the bubble. Siegel and Keshock (ref. 4) derived an 
equation in terms of the drag coefficient in the form 

F Dr=f P «°Dr ,rR iax(f) 2 < 5 > 

where C Dr was evaluated as 

_ . 45 *‘n 

/dY 

2R — (£ 

The University of Denver (ref. 20) used the same equation but evaluated the drag coeffi- 
cient as 


48 M/ 


'Dr 


2R max p *(dI) 


The larger of the two (eq. (7)) is used in this work, and the result is 


F Dr = 12 ^ R max 


(7) 


( 8 ) 
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As shown in the section Effect of Subcooling in Normal Gravity, the selection of the drag 
coefficient is unimportant in that the force was always very small. 


Dynamic Force 


The last force to be considered is termed the dynamic force F D . This force is 

y 

associated with effects on the bubble caused by the dynamics of the bubble and the liquid 
flow field surrounding the bubble. The nature of this dynamic force has not been clearly 
defined so that direct formulation is speculative. Therefore, in this work, the force is 
obtained by applying Newton's second law of motion to the generated vapor masses, or 



( 9 ) 


If it is assumed that the positive force direction is that of increasing y, as shown in fig- 
ure 1 (p. 4), the left side of equation (9) may be expanded in terms of the identified forces: 

F B + F P - F S y - F Dr + F Dy = £ ^ < 10 > 

If, at any instant, a bubble is considered to be a rigid body whose motion is described by 
the movement of its center of mass, the momentum side of the equation may be expanded 
in terms of measurable quantities such that 


'B + F P- F S - F Dr + F Dy=-f-— + V — 

B P S y Dr Dy ^ ^ dt ^ f 


V dY 


dp. 


g„ dt dt 


(ID 


When the bubble is observable on the surface, the absolute pressure changes within the 
bubble are small, so that the change of vapor density with time may be considered negli- 
gible. Therefore, when the last term on the right in equation (11) is dropped and the 
equation is solved for the dynamic force, 



Py /dV dY , y d 2 y\ 
g c ydt dt dt 2 / 


F Dr + 


( 12 ) 
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APPARATUS 


100-Foot Drop Tower 

The zero gravity experiment results were obtained in the 100-foot drop tower shown 
in figure 2. A free-fall time of 2. 3 seconds was obtained by allowing the experiment 
package to undergo an 85-foot unguided free fall. The experiment was prepared on the 
fifth floor of the tower, hoisted to the eighth floor, and suspended from the ceiling by a 
highly stressed music wire. Release of the experiment was accomplished by pressuriza- 
tion of an air cylinder that forced a knife edge into the support wire, which rested against 
an anvil. The experimental package was decelerated by allowing wooden spikes 6 feet 
long, which were mounted on the package, to embed in a box of sand 7 feet deep and 7 feet 
in diameter. 


Drag Shield 

Air resistance on the experiment package was reduced by allowing the experiment 
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(c) After test drop. 


Figure 3. - Schematic drawing showing position of experiment package and drag shield before, during, and after test drop. 









to free fall inside a protective air-drag shield as shown in figure 3 (p. 9). The drag 
shield was designed with a high ratio of weight to frontal area and a low drag coefficient, 
so that the deviation from true free fall would be minimized. As a result, the experi- 
mental package was subjected to a gravity level of less than 10“ gQ. This gravitational 
environment is termed zero gravity in this work. The sides of the drag shield were re- 
movable, so that the experiment package could be installed and removed. Prior to de- 
celeration in the sand box, the package came to rest on the bottom of the drag shield, 
which resulted in a usable zero-gravity test time of 2.25 seconds. 


Experiment Package 

The experiment package, as shown in figure 4, contained the boiling apparatus, cam- 
era and lighting equipment, power supplies, and associated controls. 

Figure 5 (p. 12) illustrates pictorially the boiling apparatus. The surface on which 
boiling was observed was a 0. 005- inch-thick, 0. 25-inch-wide Chromel strip with an ef- 
fective heating length of 0. 50 inch. The underside of the strip was mounted on an insu- 
lating material with high- temperature epoxy resin, while its ends were clamped down be- 
tween copper blocks. This assembly was mounted on a nonconducting material and en- 
closed in a glass tube. The tube was designed with a narrow neck at its top to prevent 
movement of the liquid due to surface tension effects during the transition from normal 
to zero gravity. Power to the strip was provided by a 4. 5- volt, 32-ampere-hour regu- 
lated battery supply. The surface was prepared by lapping with number 600 emery paper, 
polishing with 2/o emery polishing paper, and rinsing with alcohol. 

A second heater was mounted in the boiler and used to control the temperature of the 
bulk of the fluid. Power to this heater was supplied from an alternating-current source 
external to the experiment package and was turned off prior to a test run. 

The bulk fluid temperature was measured by a thermistor mounted on the end of a 
stainless steel probe. The position of this probe in the boiler, as shown in figure 5, was 
about 0. 5 inch above the test surface and directly over it. A thermistor was also 
mounted directly beneath the test strip in order to obtain an approximate value for the 
surface temperature. The resistance of the bulk thermistor was monitored on a digital 
ohmmeter, while the resistance of the surface thermistor was recorded with an oscillo- 
graph. 

Heat flux dissipated from the heater strip was determined by measuring the voltage 
drop across the surface with a vacuum-tube voltmeter. The resistance of the strip was 
measured accurately prior to testing with a Wheatstone bridge. 

The 16-millimeter motion picture camera, equipped with a split frame prism, pro- 
vided a filming rate of approximately 6500 pictures per second. A 35-millimeter lens, 
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Figure?. - Boiler apparatus. 


set at f4, captured a nominal area 1 inch wide by 0. 5 inch high. Timing marks were 
placed on the fil m every 1/1000 second by a timing light generator. 

Illumination was provided by four 20-watt spotlights located about the boiler in order 
to produce clear definition of the bubble outlines. The effect of heat from the lamps on 
the bulk fluid temperature over the approximately 1.25-second test duration was negli- 
gible. 
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EXPERIMENTAL PROCEDURE 

The boiler and a 2000-milliliter flask were thoroughly cleaned, flushed with distilled 
water, and polished. Distilled water was placed in the flask and boiled for 1. 5 hours to 
drive off dissolved gases. The boiler was filled with test fluid from the flask, and the 
experiment package placed in the drag shield and raised to the top of the drop tower. 
Power was supplied to the bulk heater, and the water was deaerated further for 15 to 
20 minutes. 

After the degassing processes were complete, the fluid was permitted to cool natur- 
ally, and the bulk temperature was monitored. When the bulk temperature was within 
10° F of the desired test temperature, the strip heater was actuated at a level that 
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initiated boiling. This was done to destroy any stratification that had been set up in the 
fluid during the cooling period. When the bulk temperature was within 3° to 5° F of the 
desired test value, the power to the strip heater was increased to the proper level. The 
heat flux for all experimental runs was approximately 28 900 Btu per hour per square 
foot. When the digital ohmmeter indicated that the desired bulk temperature was reached, 
the experiment package was released. The photographic equipment and lights were acti- 
vated by a system of relays approximately 1 second after release of the package in order 
to obtain data during the last 1.25 seconds of zero gravity time. Normal gravity testing 
was the same as that just described for zero gravity, except that when the desired bulk 
temperature was reached, the camera and lights were activated manually. The proce- 
dure was repeated for six bulk temperatures in each gravity field. 

The heater surface temperature measurements were made in normal gravity by 
using the aforementioned procedure. Tests were rim for a complete range of subcoolings 
with a highly polished surface and then repeated for a slightly polished one. The reason 
for testing heaters that had different surface conditions is discussed later in the report. 


DATA REDUCTION 

The bubbles recorded on the 100-foot rolls of film were viewed and measured on a 
motion analyzer that magnified the image eight times. Measurements on a bubble were 
made with the assumption that a vertical cross section was being viewed. In order to de- 
termine the position of the bubbles on the surface, however, the camera was inclined at 
approximately 5° from the horizontal so that a small error was present in the data. In- 
clination of the camera permitted selection of bubbles that were generated in the center 
of the heater and, hence, free of possible effects from the heater edges. 

In order to analyze the growth characteristics of bubbles statistically, lifetime, 
time during which a bubble is attached to a surface, and maximum radius were tabulated 
for as many as 15 bubbles at each subcooling and gravity level. Only those bubbles that 
grew undisturbed on the surface were used, so that under some conditions, fewer than 
15 were available. 

The force analysis was applied to a minimum of two bubbles at each subcooling and 
gravity level. Calculations were performed by a computer because of the large volume 
of work required. The derivatives required in equations (8) and (12) were obtained by 
the computer first fitting the second through tenth degree curves to the data points. 

Then, the curve that best represented the data was placed in the computer, which calcu- 
lated the slope of the curve at the necessary points. 

Force histories, plots of force against time, were prepared for the individual bub- 
bles so that the stimulus for removal could be studied. Force comparisons between 
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bubbles were made by a nondimens ional technique in which each positive force was di- 
vided by the total positive force and actual time was divided by lifetime. Plots of a par- 
ticular dimensionless force against dimensionless time were made for bubbles at differ- 
ent subcoolings and gravity levels. Because of the large amount of information generated 
from the procedures just described, only a representative cross section of the force data 
is presented in this work. 

Further analysis involved the selection of a critical point at which the removal forces 
imparted enough inertia to a bubble to start separation. This point was defined as that 
time during bubble collapse when the inertia of the bubble vapor, or the net force on the 
bubble, reached a maximum positive value; however, it was stipulated that the time se- 
lected should not occur during the last moments of the lifetime. These last moments 
were characterized by a rapid decrease in the bubble base diameter, so that the large 
positive inertia exhibited by some bubbles at that time was believed to be caused by, 
rather than the cause of, the impending separation. The values of the different dimen- 
sionless forces at the critical point were plotted against subcooling for normal and zero 
gravity. 


RESULTS AND DISCUSSION 
Determination of Surface Tension 

In order to calculate the surface tension force from equation (3), it was necessary 
to determine the surface tension at the heater surface. Because surface tension is a 
function of temperature, the heater surface temperature was measured. These temper- 
ature measurements were made in normal gravity, and the resulting calculated surface 
tension was used in both the normal and zero gravity force calculations. It is not to be 
inferred that the surface temperatures are not a function of gravity but rather that there 
may be some degree of error in the associated zero gravity calculations. The formula 
for the surface tension as a function of temperature was obtained from reference 25. 

The data obtained from these measurements are presented in figure 6. Each data 
point represents an average temperature taken over approximately 15 to 20 seconds. An 
average was used because the actual temperature fluctuated, as shown by a typical oscil- 
lograph trace in figure 7. Previous investigators such as Chun (ref. 11) and Han and 
Griffith (ref . 12) have indicated that boiling is dependent on the condition of the heater sur- 
face. Therefore, the heaters used for these tests were prepared so that their surface 
conditions were extremes of the surface conditions used in the photographic tests. In 
figure 6, curve I represents data for the extensively polished surface, while curve n is 
data taken for the slightly polished one. As expected, the rougher surface produced 
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Figure 6. - Effect of surface conditions and sub- 
cooling on surface temperature. Average heat 
dissipation maintained at 28 900 Btu per hour 
per square foot. 
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cooler temperatures because of the more numerous nucleation sites. The surface tem- 
peratures used to obtain the surface tension were taken from an average of the data 
points. The trend toward higher surface temperatures for higher bulk temperatures in 
the low heat flux region is also noted in reference 20. 


224, 


228 1 


Time — *- 

Figure 7. - Fluctuation of surface temperature. 


Bubble Growth Characteristics 

Investigation of the motion picture films indicated that the shape of the bubbles devi- 
ated considerably from spherical for most test conditions. This was expected from the 
results of references 18 and 19. However, an exception to this occurred at the high bulk 
temperatures, where smaller bubbles appeared to be more spherical. These results are 
due to the existence of a thermal growth layer adjacent to the heater surface and over 
which the temperature decreases from the heater surface temperature to the saturation 
temperature. The thickness of the layer can be such that all the bubbles observed, ex- 
cept the smaller ones at low subcoolings, extend out of the layer into the subcooled liq- 
uid. This causes condensation on the top surfaces of the bubbles, so that their shape 
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deviates from spherical. The fact that the smaller bubbles at low subcoolings were not 
distorted by condensation suggests that the thermal growth layer was thicker under such 
conditions. 

The bubbles observed were of two types: continuous growth, in which the volume of 
the bubbles increased, reached a maximum, and then decreased until separation from the 
surface occurred; and, oscillatory, in which the volume of the bubbles went through sev- 
eral cycles of first increasing and then decreasing volume before separation from the sur- 
face occurred. These types of bubbles have also been observed by the University of 
Denver (refs. 18 and 19). At the low bulk temperatures, all the bubbles were of a contin- 
uous growth type in both zero and normal gravities, and there was little coalescence no- 
ticeable between adjacent bubbles. As the saturation point was approached, oscillator- 
type bubbles began to appear, particularly in zero gravity. There was occasional coales- 
cence in normal gravity and extensive coalescence in zero gravity. Figure 8 presents 
interrupted photographic sequences of bubbles growing at approximately 38° and 5° F sub- 
cooling in both normal and zero gravity. The sequences are considered typical of the 
particular boiling conditions. 

As has been observed in reference 3 separation of the bubbles did not occur always 
at the heater surface; it occurred often a short distance above the surface on a shrinking 
neck, and a residue was left to form a new bubble. After leaving the surface, the bubbles 
that were generated at high subcooling at both gravity levels rapidly moved into the liquid 
bulk and condensed. In normal gravity near saturation, the bubbles left the surface more 
slowly than at high subcooling and did not condense within the field of view of the camera. 
In contrast, for low subcoolings in zero gravity, the bubbles barely separated; some re- 
turned to the surface before rising into the bulk. Several of the larger coalescing bubbles 
separated and moved so slowly into the bulk that subsequent bubbles forming on the heater 
surface were absorbed into the rising vapor mass. This latter phenomena is proposed as 
a principal heat-transfer mechanism in reduced gravity in reference 3. 

Further zero-gravity tests as close to the saturation point as possible indicated a 
large vapor mass situated above the heater strip into which the generated vapor from the 
surface was absorbed. As time progressed, the vapor mass surrounded the heater, and 
burnout seemed imminent. The primary difference, then, between saturated and slightly 
subcooled conditions in zero gravity was not separation of the bubbles but the inability of 
the vapor under saturated conditions to be removed from the vicinity of the heater sur- 
face. These results indicate the time dependence of burnout in zero gravity under satu- 
rated conditions, as proposed by Clodfelter (ref. 7). 

Bubble population decreased in both gravity levels as the bulk temperature was in- 
creased. Rough counts of the total number of bubbles generated over the test durations 
showed a moderate population decrease in normal gravity as compared with a consider- 
able decrease in zero gravity. 
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(a) Subcooling, 38.6°; gravity, normal. (b) Subcooling, 38.5°; gravity, zero. (c) Subcooling, 5.1°; gravity, normal. (d) Subcooling, 4.4 ; gravity, zero. 


Figure 8. - Bubble growing on heated surface. 
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Figure 9. - Growth of bubbles at various subcoolings. 



(a) Normal gravity. 
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<b) Zero gravity. 


Figure 10. - Rise of center of mass of bubbles at various subcoolings in norma! and zero gravity. 





Figure 11. - Effect of subcooling on bubble maximum radius and 
lifetime. 


Plots of total volume against time and center of mass against time are presented in 
figures 9 and 10, respectively. It is apparent from figure 9 that the bubble lifetimes 
were divided into two periods, growth and collapse, and from figure 10 that the transla- 
tions of the centers of mass attained greater rates early in the lifetimes when the bub- 
bles grew rapidly. 

The change in the maximum bubble radius with subcooling and gravity level is shown 
in figure 11(a). The data points represent the maximum radii of bubbles for which force 
histories and dimensionless data are presented later in this report. It can be seen that 
a range of bubble sizes was observed at any one subcooling in both gravity levels; a trend 
toward larger bubbles was evident at higher bulk temperatures. A comparison of normal 
and zero gravity data indicates that, at high subcooling, comparable values were obtained. 
As the bulk temperature was increased, the maximum radii in zero gravity became 
larger than those for normal gravity except for very low subcooling. The discrepancy at 
low subcooling was caused by the aforementioned coalescence, the result being that only 
smaller bubbles were able to grow undisturbed. The trend toward larger bubbles at 
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(b) Zero gravity. 

Figure 12. - Variation of bubble dynamic contact angle with time 
for various subcoolings in normal and zero gravity. 
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Figure 13. - Variation of bubble maximum growth rate for normal and 
zero gravity. 
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lower subcoolings substantiates the theory of a thicker thermal growth layer at increas- 
ing bulk temperatures. 

Data for bubble lifetime as a function of subcooling and gravity level exhibit much the 
same trends as do the maximum radius data, as shown in figure 11(b) (p. 19). It is not- 
able that, even though the lifetimes of only the smaller bubbles were measured at low sub- 
cooling and zero gravity, the increasing lifetime trend continues as subcooling decreases. 

The change in the dynamic contact angle with dimensionless time is shown in fig- 
ure 12 for several bubbles. The bubbles initially formed with high contact angles, ap- 
proaching 90°, decreased with increasing volume, and then increased to 90° again at sep- 
aration. The trend toward higher angles as separation is approached has been observed 
by several investigators. It was difficult to measure this angle, and it is estimated only 
to be within 10° of the true value. No clear trend regarding gravity dependence was evi- 
dent. 

The time rate of change of the volumes of bubbles, or growth rate, shows a trend to 
higher positive maximum growth rates for larger bubbles, as shown in figure 13. There- 
fore, from figure 11(a) (p. 19), a range of maximum growth rates may be expected at 
each subcooling and gravity level, and a trend toward higher maximum growth rates for 
lower subcoolings and zero gravity is predicted. 

Graham and Hsu (ref. 10) have proposed the dependence of bubble growth rate on the 
thermal layer thickness. Since previous findings in this work indicated that the thickness 
of the thermal growth layer increased as subcooling decreased and gravity level was re- 
duced, the trend in the maximum growth rates may be explained. The maximum growth 
rate was attained very early in the bubble lifetime, as has been proposed by Han and 
Griffith (ref. 12). 

A summary of the previous discussions involving the effects of subcooling and zero 
gravity on bubble characteristics is presented in table I (p. 22). 


Discussion of Forces 

Effect of subcooling in normal gravity . - Examination of the force histories in fig- 
ure 14 (pp. 23 and 24) shows that the general shapes of the pressure force and surface 
tension force curves did not change with subcooling. The magnitudes, however, did be- 
come larger as the subcooling was decreased because of the increased size of the bubbles. 
It should be noted that different scales were used on plots at different subcoolings in 
order to make the data presentation clearer. Investigation of the dimensionless plots in 
figures 15(a) and 16(a) (pp. 25 and 26, respectively) shows that the pressure force ac- 
counted for a large percentage of the total removal force at the critical point (see p. 14) 
for all subcoolings. 


21 


TABLE I. - EFFECTS OF INCREASING SUBCOOLING AND ZERO GRAVITY ON BUBBLE CHARACTERISTICS 


Characteristic 

Effect of increasing subcooling on bubble 
characteristics 

Effect of zero gravity on bubble 
characteristics 

Gravity 

Subcooltng 

Normal 

Zero 

High 

Low 

Coalesence of 

Decreased from occasional 

Decreased from extensive 

No apparent effect 

Increased from occasional 

adjacent bubbles 

to little or none 

to little or none 



to extensive 

Type of bubbles 

Continuous growth and few 

Continuous growth and 

Continuous growth 

Continuous growth and few 


oscillators to just con- 

many oscillators to 
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Figure 16. - Effect of pressure and buoyancy forces on bubble separation at various sub- 
coolings in normal gravity. 

The effect of subcooling on the buoyancy force could have been predicted from the 
maximum radius curves in figure 11 (p. 19). The dimensionless plots in figures 15(b) 
and 16(b) indicate the greater importance of buoyancy as the subcooling was decreased; 
at the critical point, buoyancy accounted for a major percentage of the total removal 
force for low subcooling as compared with a small percentage for high subcooling. 

The drag force is not shown on the force histories because it was very small; this 
force accounts for a maximum of about 4 percent of the total retention force. The exist- 
ence of this force as defined herein may be questioned if a source flow model for the fluid 
flow about the bubble, as has been proposed by Zuber (ref. 24), is assumed. For such a 
model, the liquid around the bubble flows radially away from the nucleation site. 

The role that the dynamic force played in the removal of bubbles was highly depend- 
ent on subcooling, as shown in figure 17. At the critical point, it diminished from a 
major percentage of the removal force at high subcooling to insignificance and retarda- 
tion at low subcooling. This retentive dynamic effect may account for differences that 
have been found between theoretical and experimental results regarding bubble size at 
departure. While investigating at approximately the saturation point with water, Siegel 
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(a) Bubble dynamics. 



(b) Bubble separation. 

Figure 17. - Effect of dynamic force on bubble dynamics and bubble separation at various subcoolings in normal 
gravity. 

and Keshock (ref. 3) found discrepancies between their experimental results and those 
predicted by Fritz's theory (ref. 26), an equation derived from the equivalence of surface 
tension and buoyancy. They determined that the experimental diameters were as much 
as 20 percent greater than Fritz's diameters, which indicates that something retained the 
bubbles for additional growth. Another equation for bubble size at departure, derived by 
Staniszewski (ref. 27), included a retentive dynamic term. 

In summary, in normal gravity bubbles generated at high subcooling were removed 
principally by the pressure and dynamic forces with buoyancy playing a minor role , while 
at low subcooling the stimulus for bubble removal was the pressure and buoyant forces 
with the dynamic force acting retentively. The drag force was negligible; surface ten- 
sion accounted for the retention force for all subcoolings. 

Effect of gravity level at various subcoolings . - The zero gravity force histories in 
figure 18 (pp. 28 and 29) indicate that the surface tension and pressure forces acted much 
the same as in normal gravity. Investigation of the dimensionless plots for the pressure 
force in zero gravity in figures 19(a) and 20(a) (pp. 30 and 31, respectively) shows that 
this force accounted for a large percentage of the total removal force at the critical point 
for all subcoolings. Comparison of the zero gravity curve with the normal gravity curve 
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Figure 19. - Effect of pressure and dynamic forces on bubble dynamics at various subcoolings and zero gravity. 
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(a) Pressure force. 



(b) Dynamic force. 


Figure 20. - Effect of pressure and dynamic forces on bubble separation at various subcoolings in nor- 
mal and zero gravity. 


in figure 20(a) indicates that the pressure force was not affected by the reduction in 
gravity. 

The dimensionless plots in figures 19(b) and 20(b) and the zero gravity force his- 
tories show the important role that the dynamic force assumed at the critical point for 
all subcoolings. At high subcooling, the role that the dynamic force assumed in bubble 
removal was comparable at both gravity levels, whereas at low subcooling it was not. It 
is reasonable to assume that, in order to effect bubble removal at low subcooling, the 
pressure and/or dynamic force assumed a more important role in zero gravity to make 
up for the absence of the buoyancy force. 

In zero gravity, therefore, the pressure and dynamic forces removed the bubbles at 
all subcoolings, and surface tension accounted for the retentive force. 


Dynamic Force Discussion 

In an analysis of the motion of a bubble generated on a heated surface, it is 
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necessary to recognize the existence of a force applied to the vapor mass and caused by 
the effect of the surrounding liquid set in motion by the growing and collapsing bubble . A 
search of the literature did not yield a formulation which described the phenomena accu- 
rately. In this report, Newton's second law of motion was applied to the bubble, and the 
unknown force, termed the dynamic force, was determined from the measured vapor in- 
ertia terms and the other identified forces (see eq. (12), p. 7). In this approach to the 
problem, the dynamic force term implicitly contained all the terms necessary to satisfy 
the equation of motion. Although the method of formulation yielded the magnitude of the 
dynamic force, it did not provide explicit information regarding the cause or causes of 
the force . It should also be noted that the applicability of this approach during the early 
portions of the bubble lifetime may be questioned. At this time, the growth rate of the 
bubble was very high, and the motion of its center of gravity away from the surface may 
have been dominated by the growth of the bubble in the thermal layer rather than the re- 
sult of positive applied forces. 

In considering the causes of the dynamic force, it is believed that the bulk of this 
force is associated with the momentum of the liquid surrounding the vapor; thus, a dis- 
cussion of the liquid flow fields is necessary. A flow similar to a source type, as pro- 
posed by Zuber (ref. 23), occurred for the 
growth portion of the bubble lifetime. This 
is illustrated by the bubble profiles in fig- 
ure 21(a), where it can be seen that the 
liquid- vapor interface continually moved 
away from the nucleation site. During col- 
lapse of the bubble, a source type of flow 
existed near the top of the bubble, where 
the liquid- vapor interface was still moving 
away from the nucleation site. Near the 
heater surface, however, because the bub- 
ble base was decreasing, a flow toward the 
nucleation site (a sink flow) was occurring, 
as shown in figure 21(b). 

The mechanism by which the momen- 
tum of the liquid affects the bubble has not 
been determined exactly. A proposal for 
this mechanism, made by Han and Griffith 
(ref. 12), suggests that the momentum im- 
parted to the liquid ahead of the bubble 
tends to draw the bubble off the surface if 
the growth of the bubbles decelerates 





(a) Bubble growth. 



^^/////////^^ 

(b) Bubble collapse. 

Figure 21. - Profiles of bubble at 22. 9° subcooling and zero gravity 
illustrating liquid flow patterns. 
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rapidly enough. During bubble growth, therefore, the liquid source flow and the evidence 
that the growth rate decreased rapidly provide support for this theory. Han and Griffith’s 
proposed mechanism also partially applies during collapse because of the source flow oc- 
curring at the top of the bubble. However, the sink flow near the heater surface gives 
rise to a different mechanism, one that is caused by dynamic pressure on the bottom por- 
tion of the liquid- vapor interface. 


SUMMARY OF RESULTS 

An analysis of the dynamics of a bubble generated on a flat surface resulted in a new 
formulation for the pressure, buoyancy, and dynamic forces. Application of the analysis 
to experimental data taken from high-speed motion pictures of boiling water at low heat 
flux and various subcoolings in normal and zero gravity indicated the following results: 

1. The newly defined pressure force was of major significance and accounted for a 
large percentage of the total removal force at the critical point for all subcoolings and 
gravity levels. 

2. In the discrete bubble region at high subcooling, the pressure and dynamic forces 
as defined herein were the primary stimuli for removal in both normal and zero gravi- 
ties; this accounts for the gravity- independent nature of boiling phenomena exhibited at 
such conditions. 

3. At low subcooling, the importance of the dynamic force in bubble removal in- 
creased in zero gravity as compared with normal gravity, which makes up for the ab- 
sence of the buoyancy force and, hence, aids bubble removal. 

4. The drag force is of minor importance in the retardation of bubbles growing on a 
surface. 

Observation of boiling at the various test conditions and measurements of bubble 
characteristics resulted in the following conclusions: 

1. The thermal growth layer thickness appeared to be a function of subcooling and 
gravity level, except at high subcooling, where the layer seemed to be independent of 
gravity level. 

2. Continuous -growth bubbles were generated at all subcoolings and gravity levels, 
and some oscillatory bubbles appeared as subcooling was decreased in normal and zero 
gravities. 

3. Bubble maximum radius and lifetime were functions of subcooling and gravity 
level, except at high subcooling, where they were independent of gravity level. 

4 . The bubble dynamic contact angle approached 90° at bubble inception and separa- 
tion and did not appear to be affected by subcooling and gravity level. 



5. The growth rate of a bubble reached a positive maximum early in its lifetime and 
appeared to be a function of the thermal growth layer thickness. 

6. In zero gravity when the liquid temperature was close to the saturation point, the 
generated vapor lingered in the vicinity of the heated surface. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 8, 1966, 

124-09-03-01-22. 
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APPENDIX A 


BUOYANCY FORCE CALCULATIONS 


The net upward differential pressure 
areas, as shown in figure 22, is 



force acting on the differential bubble surface 
dFg = Pj cos dAj^ - P 2 cos 0 2 dA 2 (^1) 


cos 6 ^ dAj = cos 6 2 dAg = 2?:r dr (A2) 


p i - p 2 = Pe 7- (?2 - *1) (A3 > 

S C 

Substituting equations (A2) and (A3) in (Al) 
yields 

dF B = Pf — (y 2 - yi ) 2vr dr ( A4 ) 

g c 

Integrating over the surface area that is not 
directly over the base results in 


F B = f ( y 2 " y l) 2lir dr (A5) 

Sc *A-A b ’ 

However, 

/ (y 2 - yx)2irr dr 

A-Ag 


is the total volvune of the bubble minus the volume directly over the base. Therefore, 


F B'P«f-( V - V b) 

g c 


However, 


and 
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APPENDIX B 


PRESSURE FORCE CALCULATIONS 


The net upward pressure force acting on a differential bubble surface area directly 
over the base (fig. 22, p. 35) is 

dFp = (P. - P Q j cos a dA (Bl) 

Assuming that the top of the bubble is a segment of a sphere with radius of curva- 
ture T, 



(B2) 


(ref. 28) and 

cos a dA = 2irr dr 

Substituting (B2) and (B3) in (Bl) yields 



(B3) 


(B4) 


Integrating gives 

„ 2 "W t b 

*p 

T 

which can be expressed in terms of the base diameter as 



(B5) 


(B6) 
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A motion-picture film supplement C-246 is available on loan. Requests will be filled 
in the order received. You will be notified of the approximate date scheduled. 

The film (16 mm, 11 minutes, color, sound) shows the effect of subcooling and grav- 
ity level on the characteristics of water vapor generated on a flat horizontal surface. 

Film supplement C-246 is available on request to 

Chief, Technical Information Division (5-5) 

National Aeronautics and Space Administration 
Lewis Research Center 
21000 Brookpark Road 
Cleveland, Ohio, 44135 


i _CUT__ 

[ Date 

! Please send, on loan, copy of film supplement C-246 to 
!TND- 

iName of organization 
: Street number 

| City and State Zip code 

! Attention: Mr. 

! Title 


"The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.” 

— National Aeronautics and Space Act of 1958 
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